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Abstract 

We present parametrizations for the proton structure function F2 in the next to leading 
order in perturbative QCD. The calculations show that the dominant term to F2{x,Q'^) 
should grow as x~'^^ for small x values, with the exponent A5 being essentially independent 
of Q^. Comparisons with the most recent HI and ZEUS data confirm the value Ag ~ 0.35 
obtained previously from fits to low energy data. 
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1 Introduction 



One of the most interesting results obtained at HERA so far, has to do with the dramatic 
increase of the proton structure function F2 at low Bjorken x values |]^ and |^. The 
precision of these measurements do allow the extraction of the gluon density in the proton 
down to X ~ 10"^ and ||] . 

Although the general framework to discuss deep inelastic scattering is established and well 
known since two decades, the interpretation of these data is subject to some controversy. 
Two lines of thought are generally followed to describe the data. 

• On the one hand, there are those who advocate that this dramatic increase of the 
proton structure function can be obtained from singular H, ^ (non-singular 0) 
parton densities at moderate Qq ^ A GeV^ (resp. at very small Qq ~ 0.5 GeV^), 
which are then evolved using the well known DGLAP equations Q. This procedure 
describes the experimental data very well at the cost of having approximately twenty 
parameters which enter into the parametrization of valence and sea quark, as well 
as gluon densities at an input Qq value, in addition to assumptions about their 
functional forms. 

• On the other hand, there are those who argue that since at very low x values, 
the boson gluon fusion mechanism is the dominant source of leading order (LO) 
corrections to the Born level cross sections and the kernel Pg — > gg is singular, one 
expects that when including higher order QCD corrections, one would have to sum 
multigluon exchange ladders. Depending on the approximation used to perform this 
sum, one encounters deviations from the DGLAP linear evolution equations. In fact, 
the proponents of this approach claim that at fixed the x-dependence is most 
generally given by the BFKL equations |^, which predict the parton densities to 
behave like x~'^ with 

vr 

for ~ 20 GeV^. However, theory does not tell us where in the {Q'^,x) plane 
the transition region lies beyond which the expansion in terms of log ^ is impor- 
tant. The BFKL evolution equations are not yet numerically implemented in the 
parametrizations discussed above. 

We think that in order to clarify the issue one should try 

• to search for specific fingerprints of the BKFL equations, as discussed by Mueller 



[10 1 and others |11]. In particular since the above prediction runs contrary to the 



trend observed in the data |12| 

to make detailed comparisons between analytic NLO QCD predictions and experi- 
mental data to see if one can isolate regions of phase space where discrepancies might 
appear. In this context, we would like to remind the reader that since the behaviour 
at small x of the proton structure function is connected with the singularities of 
the operator product expansion matrix elements, one has two specific predictions, 
depending on whether these singularities lie to the left or to the right of those of the 
anomalous dimension matrix [12|, such that either 



— the proton structure function should grow faster than a log but slower than a 



power in 2; [13|, i.e 



F2(x,Q2) - Co 



33 - 2n 



/ 



1 
4 

exp , 



5767r2[logx|log[a,(Q2)/«^(Q2)]J ^^(33-2n/) [ ^ a.(Q2) 
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leading to a double asymptotic behaviour as discussed by Ball and Forte 10 

— or else the proton structure function should behave as a power in x i.e. x"'*'^ 
with A5 a independent constant except for possible variations due to crossing 
flavour thresholds which could depend on [jl^, and |16]. This behaviour 
is the same as that calculated by Witten for 7*7 scattering [p^]. 

The purpose this paper is to test this second set of predictions, which date as far 
back as 1980 pi and dwell upon the possibility that the cross sections for off- 



shell particles grow as a power of centre of mass energy PU^ , see also |26]. This is 
particularly interesting, since it has been shown that although the double asymptotic 
behaviour is a dominant feature of the data, there appear non-negligible scaling 



violation effects [22|. 
Our aim is twofold 

* to present NLO parametrizations for F2 and make extensive comparisons with the 
most recent data and [p^]. 

* to extract the gluon density in the HERA kinematic regime and to present predic- 
tions for R{x, Q"^) which will be measured soon at HERA. 



2 LO and NLO predictions 



Perturbative QCD provides evolution equations for the structure functions, in such a way 

that if we know them at a given Qg, we can predict them at any other value, assuming 

that both Qg and lie in a range where perturbation theory applies. 

The DGLAP equations represent one of the forms in which the evolution in can 

be expressed, whereas the most direct result of the operator product expansion (OPE) 

approach is expressed in terms of moments of structure functions. Both approaches are 

equivalent. 

Once we know the evolution equations, the goal will be to find the functional form F(x,Q'^) 
for the structure functions such that used at an input Qq, the resulting 'evolved' function 
would continue to be the same F{x, Q^) simply calculated at the new value. No such 
functional form has been found so far, thus any simple analytical form chosen as input is 
not invariant as a function of Q^, in contradiction with the fact that any other value Qf 
could have been chosen as the starting point for the evolution. 

However, one can find functional forms compatible with QCD in definite x regions; so 
we can give simple functional forms compatible with the DGLAP evolution equations, in 
particular the behaviour at the end points x = 0,1 and certain sum rules. As a result, we 
know exact solutions to the DGLAP evolution equations, albeit only locally and not for 
the whole x range. These results were derived in [19] and |24] at the leading (LO) and 
next to leading (NLO) order. 

Let us consider F2{x, Q^) for DIS e — p scattering. It can be written as 

F2{x,Q^) = Fs{x,Q^) + Fns{x,Q^) (3) 

where Fs{x,Q'^) {F^si^^Q'^)) refer to the singlet (resp. non-singlet) pieces, whose evo- 
lution equations are different. Indeed the singlet part and the gluon momentum density 
Fg{x,Q'^) evolve together, whereas the non-singlet part is decoupled from gluons and 
evolves independently. 
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2.1 LO predictions 

To LO, the evolution of the moments takes a very simple form 



IJ'Ns{n,Q^) = 



2^ d(n) 



for the non-singlet piece and, for the singlet, 

'asiQo 



■ fJ-Nsin, Ql) 



(4) 



2\ D(n) 



jl{n,Q 



il{n, Ql) 



(5) 



where /2(n, Q^) is the two-component vector 



( /^5(n,Q2) \ 



defined as usual 



l_i{n,Q^)= [\xx''-^F{x,Q^) 
Jo 



(6) 



(7) 



and d(n) and D(n) are proportional to the anomalous and anomalous dimension matrix, 
whose exact expressions are 



dfn) 



16 



+ - - Si[n, 



33 -2n/ \2n{n + l) 4 



(8) 



and 
D(n) 

where 



16 



33 -2n/ 



2^ + -4-Si{n) 

n^+n+2 
2n(n^-l) 
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4n(n-l) + 4(n+l)(n+2) 



3n/ n^+n+2 
8 n(n+l)(n+2) 
9 33-2ny _ 9Si(n) 

16 1 



5i(n) = n^ 



(9) 
(10) 



and is the number of flavours. 

Prom the fact that Qq and are arbitrary values in the range of applicability of pertur- 
bation theory, it follows that the moments must be of the form 



2 r 2 1 ~'^(") 

ms{n,Q ) = asiQ ) -hNsin) 



and 



-D(n) 



6(n) 



with b^sin) and &(n) being independent of the squared momentum transfer. 
If we try a Regge inspired functional form 

FNsix,Q'') = BMsiO") ■ x^^^^'^"^ 



and 



-As(Q") 



(11) 
(12) 

(13) 
(14) 
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for the structure functions near the end point x = 0, the dependence of B{Q'^), X{Q'^) 
upon is related to the moments in Eqs. 11 and 12 for the value of n at which they 
diverge. 

The results which can be found in |24] and |1£] are such that Xns and A should be 
independent. Furthermore one must have Aat^ < 1 and > 0, and 



and 



Bns{Q^) = Bns- as{Q' 



Bs{Q'') = Bs- asiQ^ 



-d(l-Ajvs) 



(15) 



(16) 



In addition, the gluon structure function ought to be proportional to the singlet piece of 
the structure function i.e. 

Fg{x,Q^) = Bg{Q^)-x-^s (17) 



with 



d+{l + Xs) - Dn{l + Xs) 



BsiQ^) = Bgs ■ BsiQ^ 



(18) 



^12(1 + As) 

where by d^{n) we denote the largest eigenvalue of the anomalous dimension matrix D(n). 




Figure 1: The dependence on A of Bgs^vg and vs 



2.2 NLO predictions 

The extension of these results to the NLO is tedious, the result being of the form 



Bns{Q^) = Bns ■ [«s(Q') 
BsiQ^) = Bs ■ lasiQ^] 



• |H ^ vns{^-Xns)\ (19) 

■{1 + ^^-vsil + Xs)} (20) 



and 



-d+{l+Xs) 



{Bgs + 



vg{1 + Xs)} 



(21) 



with VNsi^ ~ ^Ns), vs{l + A5) and + ^s) being known functions of the exponents 
\ns and \s given in H Q. 

The last two of these functions have been fitted to simple rational expressions and the 
results for n/ = 3,4 are presented in Fig. 1 together with similar fits to Bgs which is 
defined in Eq. 18 as the proportionality factor to LO between Bg{Q'^) and Bs{Q^)- 
An important comment to make is that, NLO corrections to the singlet and gluon structure 
functions are large, in contrast to the situation for the non-singlet case. To be rigorous 
one would have to replace Eqs. 20 and 21 by the more precise exponential forms that 
follow from the evolution equation, and of which Eqs. 19-21 are an approximation. 
Note that the gluon structure function, at the end point x = 0, is completely determined 
by the quark singlet structure function which in turn is the dominant piece of F2{x^Q'^). 
In fact, 

Bgs + j^-vg{1 + Xs) 



Fg{x,Q' 



l + f--vs{l + Xs) 



Fs{x,Q' 



Since Fs{x, Q"^) is defined in terms of quark densities in the proton as 
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(23) 



where the sum runs over quark flavours with charge given by e^. The gluon density is 
then given by 

xG{x, Q2) = — L_ . Fg{x, g2) (24) 
< > 



Let us now turn to a discussion of the limit 2; ^ 1. We have considered the ansatz 



F(x,Q2 



A{Q'')-{l-x 



(25) 



The dependence of the functions A{Q'^) and fiCp) is now related to the evolution of the 
moments in the limit n — > 00. One finds 



with 



and 



Ans{Q^) = Ans- [as(Q' 

^s(Q') = As ■ [as{Q^) 
16 



-do 



HQ' 



33 - 2n 



r(i + up) 

r(i + 1^0) 
■r(i + KQ2)) 

log as{Q^) 



do 



16 .3_ 
33 -2n/ U 



(26) 
(27) 

(28) 

(29) 



with Euler's constant. 

The gluon momentum density is fully determined, also in this limit, by the singlet structure 
function: 

Fg{x,Q^) 



2 2 (1 - x)''('3^)+^ 



(30) 



^Note that there are a few typographical errors in this reference, they will be corrected in a forthcoming 
publication by K. Adel and F.J. Yndurain 
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Finally from the fact that d{l) = and d-^-{2) = 0, the following sum rules can be derived 



• Fns{x, Q^)dx = (l - 3 < ef >) {1 + ^as(Q')} 



(31) 



with k a known and very small constant |24| and 

j\Fs{x, Q^) + Fcix, Q^)}dx = < > .{1 - ^^siQ^)} (32) 

3 Parametrizations for the proton structure func- 
tion 

In this section, we give approximate parametrizations compatible with the exact conditions 
discussed above. 



3.1 LO parametrizations 

We consider the following ansatz for 



Fsix,Q^) = {BsiQ^)-x~'^^ + CsmUl - xyC^'^ 



where 



Bs{Q^) = Bs- as{Q 



-d+(l+\s) 



(33) 



(34) 



and v{Q'^) given by Eq. 28. 

The singlet part is the most important contribution at small x, so that the predicted 
behaviour near x = has to be implemented through the term BsiQ"^) • x~^. Note from 
Eq. 33 that it is divergent when x — > and that the coefficient BsiQ"^) grows rapidly 
with Q^, becoming asymptotically dominant. This part is responsible for the increase of 
the structure function at small x as a function of as observed experimentally. 
On the other hand, in the limit —^ , Bs{Q'^) becomes negligible, in fact, it vanishes for 
= 0. At small we expect soft Pomeron-like contributions to remain. This term we 
parametrise phenomenologically with the help of the second piece in Eq. 33 proportional 
to CsiQ"^), whose relative weight will decrease as a function of increasing 

but remain 

important at low and intermediate values. 

We would like to remark that the dependence of both functions Bs{Q'^) and Cs{Q'^) can 
be determined by implementing the limits at x = 0, 1. Thus, 



Bs- a.(g' 



-d+{l+\s) 



+ As- as{Q 



-do r(i + i/o) 



For the non-singlet part we could try a parametrization of the form 



FNsix,Q^) = {BNsiQ^ 



+ CNsiQ'') • x} • (1 - x) 



HQ'') 



but since its contributions turns out to be small, we will use simply 



Fns{x, Q2) = Bns{Q^) ■ x-^^^ ■ (1 - x) 



with 



Bns{Q^) = Bns as{Q 



-do r(l + i/o) 

r(i + KQ')) 



(35) 

(36) 

(37) 
(38) 



7 



In Eq. 36, X^s is related to the intercept of the leading Regge trajectory (p) contri- 
bution to this piece, Xns = 0.45. In one wishes to restrict the number of free parameters, 
one can use Eq. 31 to obtain 

Bns{Q ) = (1 - 3. < e, >j . f(X^^^)r(rTKQ^ ^''^ 

For the gluon momentum density, the prediction following from our analysis will be 

Q2) = {Bg(Q') ■ x-^^ + Cg(Q')} • (1 - (40) 

3.2 Parametrizations at the NLO 

At the next to leading order, the only significant modification to the expressions given 
above refer to an extra term of the form 1 + ^^/o5^(l — x) which will be significant only 
at large x. Therefore we have for the dominant singlet part 

Fs{x, Q2) = {5^(q2) . ^^Xs + Cs{Q''m - xr^'^'^ ■ {1 + ^lo9\l - x)} (41) 
with v'i{Q^) being given by the following expression 

^i(Q') = KQ') - asi^H'^iQ^) + 1) + «i} (42) 

The coefficient BsiQ'^) is given by Eq. 19, and Cs{Q'^) is such that Eq. 35 is satisfied i.e. 
including NLO terms 



-d+{l+\s) 



-do e^'("=)"'' • r(i + 



Cs{Q') = -Bs- [as{Q')\ + As ■ [as{Q')\ ■ ^(-^^^(^2)) (43) 

with 

g{as) = ao + ai • ^(1 + i^{Q^)) + + i^iQ^) - + z^(Q'))} (44) 

Here ip denotes the logarithmic derivative of the T function, oq and ai are constants 
essentially independent of the number of flavours. Numerically oq = —1.18 and oi = 0.66. 
For the non-singlet part we take 

Fns{x, Q') = BNsiQ') ■ x-''^' ■ (1 - x)^(«') • {1 + ^log^l - x)} (45) 

OTT 

with Bf^siQ"^) fixed in such a way so as to satisfy the sum rule given in Eq. 31, whose 
NLO corrections are negligible. These parametrizations are similar to those used by the 



authors of |25| to fit fixed targed lepton nucleon scattering data. At this point we would 
like to remark that these predictions are valid for a range in where the number of 
flavours is fixed. 



4 Comparison with experimental data 

We can think of two possibilities to compare experimental data to theoretical predictions. 
The simplest is illustrated in our fits of the '94 ZEUS shifted vertex data [pll] to the 



expression given by Eqs. 14 and 20 in NLO. This has also been tried in |12]. In order to 
be in the kinematical region where this term dominates the singlet structure function, we 
restricted ourselves to values above 3 GeV^ and Bjorken x below 0.01. The QCD scale 
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parameter A^j^. was fixed in a way be discussed below. The reason for this has to do with 
the fact that the normalization factor Bs and the strong coupling constant are tightly 
correlated. Therefore we prefer to fix the coupling constant to values measured elsewhere. 
The number of excited flavours was assumed to be constant and the quality of the fits 
are similar irrespective of whether we take n/ = 3ornj = 4. Disregarding the overall 
normalization factor, Bs, the relevant fitted parameter A5 turns out to be 0.32 ±0.01 for 
four excited flavours. The is 25 for 33 experimental points. The agreement between 
data and NLO QCD predictions is quite good as illustrated in Fig. 2, even down to small 
values where the applicability of perturbation theory could be questioned. The solid 
lines in Fig. 2 represent the results of the fits. 

The gluon density extracted from these fits also agrees well with that obtained by ZEUS 




Figure 2: The 7* — p cross-sections as a function of W 

using the Ellis-Kunszt-Levin method [^], which relates the gluon density to the loga- 
rithmic derivative of the proton structure function, following a numerical approximation 



originally derived by [28|. Our results are shown in Fig. 3. Our smaller error band has 
to do with the fact that the gluon density determination within our formalism, is subject 
essentially only to the statistical errors of the proton structure function itself. 
The criticism to the procedure developed so far is clear from Fig. 2. A fast rise of the cross 
section for virtual 7-p scattering is observed at large W'^ values, but this rise is sitting 
on top of a non-negligible plateau which moreover exhibits clear leading twist behaviour. 
Questions like, what is the dependence of A5 upon variations in the limits used to define 
the kinematical region over which the fits are done, or whether one should subtract from 
the low X region a contribution smoothly coming from the large x (i.e. small W'^) domain, 
have to be answered in a quantitative way. 

In order to do this, we consider a different approach in which we have performed fits to 
the data in the whole {x,Q'^) region covered by the experiments. We would like to recall 
that four parameters are involved in our fits, if we consider a region in with a fixed 
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number of flavours, namely a coefficient Bs to give the normalization of the singlet piece 
to F2, a coefficient As which serves to define the subleading contribution to the singlet 
piece, Xs which defines the growth rate for small x and vq which fixes the behaviour of 
the structure function at large x values. The strong coupling constant is fixed through 
^MS ~ MeV for four flavours as given in |2^. The dependence of the QCD scale 
parameter on the number of flavours is taken as in [^]. Although reasonable fits can be 
obtained with nj = 4 independent of Q^, we find that the quality of the fits improve by 
considering n/ = 3 below = 10 GeV'^, n/ = 4 for between 10 GeV'^ and 100 GeV'^ 
and Uf = 5 above. In principle, Xs could be different for different number of excited quark 
flavours. Therefore, we have tried two sets of fits, one with Xs fixed for different quark 
flavours and a second one allowing different Xs values as a function oinj. The second set 
of fits yield considerably improved values, and these are shown in the figures. We are 




Figure 3: The gluon density determined from shifted vertex ZEUS data, dotted band, compared 
with previous determinations. 

aware that this is an approximation: for instance in the model of GRS the charm quark 
contribution rises smoothly with increasing Q^. Once F^^"-"^"^ is measured, this could be 
subtracted from the experimental F2 values and we could keep n/ fixed to = 3 in the 
complete range. The effect of the bottom quark excitation, due to its smaller charge, 
is smaller. These are the main sources of uncertainty in our determination of the gluon 
density and of R{x,Q'^) which we cannot take into account in a model independent way 
right now. 

The HI data |^] has been fitted over the entire x region and for larger than 3.5 GeV^. 
In Fig. 4 we show the results of the fits. We also show the extrapolation from our fits to 

down to 1.5 GeV^ where the validity of the perturbative expansion is debatable. The 
gluon density extracted from these fits is presented in Fig. 5 for a restricted range 
between 3.5 and 250 GeV^. One should note that in the first two bins, the predictions 
tend to be above the data systematically. There could be several reasons for this: miss- 
ing next to NLO corrections, only two flavours in the proton are excited, or more likely 
non-perturbative effects, such as higher twist, begin to emerge. 

We would like to point out that at very high Q2 values, i.e. > 800 GeV^, the phe- 
nomenological term Cs{Q'^) which we attribute to the soft pomeron contribution, tends 
to become small and negative, an effect due to the oversimplified parametrization over 
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Bs 




As 




/n.o.p. 


HI 


(5.4 ±0.3) • 10"^ 


0.36 ± 0.01 


0.46 ±0.05 


1.3 ±0.1 


162/181 


ZEUS 


(4.6 ±0.3) • 10-4 


0.35 ± 0.01 


0.53 ±0.04 


1.2 ±0.1 


240/175 



Table 1: Summary of the fitted parameters 





Bs 


As(Q2 < 10) 


As(10 < Q2 < 100) 


As(Q2 > 100) 


As 




X^/n.o.p. 


HI 


(4.2 ±0.3) ■ 10^-* 


0.33 ±0.01 


0.33 ±0.01 


0.36 ±0.01 


0.40 ± 0.04 


1.1 ± 0.1 


112/181 


ZEUS 


(3.5 ±0.3) • lO""* 


0.32 ±0.01 


0.33 ±0.01 


0.35 ±0.01 


0.42 ± 0.03 


1.0 ±0.1 


202/175 



Table 2: Summary of the fitted parameters 



the whole (x, Q^) range rather than to a genuine physical effect. 

The ZEUS data [^] has been fitted over the entire x range for larger than 6 GeV^. 
The results of the fits are presented in Fig. 6 up to = 650 GeV^, and the corresponding 
gluon densities up to = 200 GeV'^ in Fig. 7. 



The values of the parameters obtained from fitting the HI |22| and ZEUS data |23] 
are summarized in Tables 1 and 2. They are reasonably similar. The for the fits of the 
ZEUS data are poorer. The main contribution to the comes from the high data. 
Any attempt to fit a smooth function to these data is bound to give poorer values. We 
have tried some variations in the fitting procedure by considering the phenomenological 
term CsiQ"^) to be proportional to x'', with e an additional free parameter. The resulting 
fitted value for e turned out to be within errors compatible with 0. 

We do not expect the parameter z^o to be well determined by the fits, since it is linked 
to the behaviour of the structure functions at large x values, while the HERA data are 
concentrated in the low x domain. A more reliable determination of i^q would require to 
incorporate the data from fixed target experiments. 

5 Predictions for R(x, Q^) 

As it is well known, in the quark parton model the Callan-Gross relation holds, namely 

F2{x,Q^) = F,{x,Q^) (46) 

in such a way that the ratio 

, _ F2{x,Q^)-Fi{x,Q') _ FLix,Q') 

Fi(x,Q2) - F,{x,Q') ^^^> 
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vanishes. QCD corrections induce violations of the Callan-Gross rule, leading to non-zero 
values for R{x,Q'^). Defining the non-singlet and singlet contributions to Fl as 



and 



with 



dry r 

Fl,ns{^,Q') = ^ dy-F2,Ns{y,Q') (48) 

FlA^, Q2) = ^{ dy^F^Ay, Q') + C dy^{l - -)FG{y, Q^)} (49) 
37r Jx y-^ Jx y 

5f = l-nf (50) 



and using the parametrizations discussed in previous sections for Fs and Fq, one obtains 

This expression, originally derived in |31] should be valid for large values, where the 
Pomeron like contribution is negligible. At low and intermediate values, one has to 
consider the full parametrization for both Fs and Fq^ so that one obtains 

Notice that both expressions are proportional to Og. In Fig. 9 we present the predic- 
tions from Eq. 52, lower band, derived from the fits to the ZEUS data. The predictions 
for R exhibit for a given range a smooth rise with decreasing x up to a value of ap- 
proximately 0.25. In order to indicate the importance of the Pomeron-like term, we also 
show the predictions derived from Eq. 51, upper band, which was calculated by setting 
CsiQ'^) = 0- We expect future measurements of R to lie close to the lower values. 



6 Conclusions 

We have compared recent HERA data on structure functions at low x with QCD analytical 
calculations. NLO predictions in QCD describe the rate of growth of the proton structure 
function F2 in a wide and x domain. The dominant term behaves like x~'^^ with 
Xs ~ 0.34 lb 0.03 independent of Q^. This spread takes into account possible dependences 
on the number of excited quark flavours. This is in contrast with recent results by the HI 
Collaboration which suggested that A5 grows from 0.08 at low to 0.5 at high Q'^. We 
can exclude a BKFL prediction where the exponent uj in the x behaviour is proportional 
to the strong coupling constant as{Q'^) and therefore decreases with Q^. 
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Figure 4: The HI 1994 data along with the results of the fits described in the text 
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Figure 5: The gluon density extracted from the HI 1994 data. The dotted bands indicate the 
uncertainties in the fitted parameters. 
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Figure 6: The ZEUS 1994 data up to ^ 650 GeV^ along with the results of the fits described 
in the text. 
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Figure 7: The gluon density extracted from the ZEUS 1994 data. The dotted bands indicate 
the uncertainties in the fitted parameters. 
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Figure 8: Predictions for R{x,Q'^) from fits with CsiQ"^) equal and non equal to zero, as 
described in the text. The shaded bands represent the uncertainties in the fitted parameters. 
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